Coexistence of spherical states with deformed and superdeformed bands in doubly 

magic ""^Ca; A shell model challenge. 
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Large scale shell model calculations, with dimensions reaching 10^, are carried out to describe 
the recently observed deformed (ND) and superdeformed (SD) bands based on the first and second 
excited 0^ states of ''"Ca at 3.35 MeV and 5.21 MeV respectively. A valence space comprising 
two major oscillator shells, sd and pf, can accommodate most of the relevant degrees of freedom 
of this problem. The ND band is dominated by configurations with four particles promoted to 
the p/-shell (4p-4h in short). The SD band by 8p-8h configurations. The ground state of ^''Ca is 
strongly correlated, but the closed shell still amounts to 65%. The energies of the bands are very well 
reproduced by the calculations. The out-band transitions connecting the SD band with other states 
are very small and depend on the details of the mixing among the different np-nh configurations, 
in spite of that, the calculation describes them reasonably. For the in-band transition probabilities 
along the SD band, we predict a fairly constant transition quadrupole moment Qo (1)^170 e fm^ 
up to J=10, that decreases toward the higher spins. We submit also that the J=8 states of the 
deformed and superdeformed band are maximally mixed. 

PACS numbers: 21.10.Sf, 21.60.Cs, 23.20.Lv, 27.40.+Z, 29.30.-h 
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I. INTRODUCTION 

^°Ca is a textbook example of doubly magic nucleus, 
corresponding to neutrons and protons filling the three 
lowest harmonic oscillator shells with principal quantum 
numbers p=0, 1, and 2. The N=Z=20 gap is large (about 
7 MeV) and therefore, in a naive independent particle 
description, its spectrum will consist of negative parity 
quasi degenerated states starting at an excitation energy 
commensurable with the value of the gap and positive 
parity states at about twice the value of the gap. In- 
deed, this is not the case and its first excited state is 
a 0+ at 3.35 MeV, that turns out to be the head of a 
deformed band. Excited deformed bands in spherical nu- 
clei provide a spectacular example of coexistence of very 
different structures at the same energy scale, that is a 
rather peculiar aspect of atomic nuclei. Other cases have 
been known for a long time, for instance, the four parti- 
cle four holes and eight particle eight holes states in ^^O, 
starting at 6.05 MeV and 16.75 MeV of excitation energy 
P, 01 . The theoretical descriptions based in multiple par- 
ticle hole excitations that can accommodate deformation 
started with the work of Brown and Green Q, Gerace 
and Green 0, [i| and Zuker, Buck and McGrory Q In 
^"Ca, there has been since long experimental indications 
showing that the two low- lying sequences, 0+, 2+, 4+, 
starting at 3.35 MeV and 5.21 MeV, may correspond to 
deformed or superdeformed bands (see ref. 0). How- 
ever, it is only recently that such bands, deformed and 
superdeformed, have been explored up to high spins, in 
several medium-light nuclei such as ^^Ni Q, ^^Ar [91 and 
of course ''"Ca |10|], thanks to the availability of large 
arrays of 7 detectors like Euroball and Gammasphere. 

One characteristic feature of these bands is that they 
may belong to rather well defined shell model configu- 



rations; for instance, the deformed excited band in ^^Ni 
can be associated with the configuration (Of7/2)^^ (lP3/2j 
0^5/21 lPi/2)'' while the super-deformed band in '^^Ar is 
based in the structure (sd)^^ (pf)'^- The location of the 
np-nh states in "^^Ca was studied in the Hartree-Fock 
approximation with blocked particles and Skyrme forces 
in ref. [llj . While many approaches are available for the 
microscopic description of these bands (Cranked Nilsson- 
Strutinsky [o'l, Hartree-Fock plus BCS with configura- 
tion mixing ^12J . Angular Momentum Projected Genera- 
tor Coordinate Method [H, Projected Shell Model [11], 
Cluster models [H, [11] etc.) the interacting shell model 
is, when affordable, a prime choice. The mean field de- 
scription oi N = Z nuclei, has problems related to the 
correct treatment of the proton-neutron pairing in its 
isovector and isoscalar channels, as well as the proper 
restoration of the symmetries broken at the mean field 
level to take into account the pairing (particle number) 
and quadrupole (angular momentum) correlations. In 
addition, the triaxial degrees of freedom are seldom con- 
sidered. On the shell model side, the problems come from 
the size of the valence spaces needed to accommodate the 
np-nh configurations. 



In a recent article [17[, it was shown that large scale 
shell model calculations within a valence space consisting 
on the sd and pf shells can describe the coexistence of 
the spherical ground state and the superdeformed band 
of ■^^Ar, both the very existence of the spherical and the 
superdeformed states forming a band and the configu- 
ration mixing that permits transitions connecting them. 
In the same spirit, in the present work we focus on the 
case of ^'^Ca, which is particularly challenging since it 
presents two rotational bands of different kind above the 
spherical ground state. The first rotational band, start- 
ing at 3.35 MeV of excitation energy, is generated by an 
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intrinsic deformed state, whereas the other one, starting 
at 5.21 MeV above the ground state, stems from an in- 
trinsic super-deformed state, i.e. bearing an axis ratio 
close to 2:1. 

In earlier works we have surmised that the superde- 
formed band corresponds to a structure (sd)^^ {pf)^, 
while the normal deformed band presents a (sd)^" (pf)'^ 
structure [l^, U^. That is, the structure of the bands 
corresponds to the promotion of four and eight particles 
across the Fermi level, from the sd to the pf shell, that 
we call 4p-4h and 8p-8h configurations. However, this de- 
scription is too crude, because the physical states contain 
components of different np-nh rank. The mixing should 
allow for transitions connecting the superdeformed and 
normal deformed bands as well as both bands with the 
spherical ground state, while remaining gentle enough not 
to jeopardize the very existence of the bands. This is the 
aim of this work; to explain the coexistence of spherical, 
deformed and superdeformed structures at low energy in 
""^Ca and to understand how do their dominant configura- 
tions mix among themselves while preserving their iden- 
tity. 



II. SPHERICAL MEAN FIELD vs 
CORRELATIONS 

The intriguing question is how comes that these many 
particle-hole configurations appear so low in excitation 
energy. Before answering it with the actual calcula- 
tion, we want to advance a more qualitative discussion 
based in the competition between spherical mean field 
and quadrupole and pairing correlations. 

In a very recent paper 20] Rowe et al. have ex- 
amined this issue in ^^O, with a model consisting of a 
spherical harmonic oscillator mean field plus an Elliott's 
quadrupole-quadrupole interaction. In their model, the 
unperturbed excitation energy of an np-nh excitation 
between contiguous shells is nhuj . With the stan- 
dard prescription for hui the value of the ''°Ca gap is 
largely overestimated (12 MeV instead of the experimen- 
tal 7 MeV). The quadrupole correlation energy is taken as 
proportional to the Casimir operator of the lowest grade 
SU(3) irrep corresponding to the np-nh configuration, 
C's(7(3): with the coupling constant for the quadrupole- 
quadrupole interaction equal to hui /2No, (Ng being the 
total number of oscillator quanta in the ground state of 
the nucleus of interest). The excitation energy of the 
np-nh band-head can be written as: 

AE = (^(N(np - nh) - No) - (1) 

For the deformed band -n=4, A=12, ^==8- AE=1.06fiw 
and for the superdeformed band -n=8, A=20, ^=12- 
AE=0.66huj (A and fi being Elliott's labels of the irreps 
of SU(3)) This model is very schematic, and, as we shall 
discuss later, overestimates both the energy differences 



between the closed shell and the np-nh configurations 
at the mean field level, and the correlation energies of 
the later. Thus, it is not surprising that it cannot repro- 
duce the experimental numbers. But it contains the right 
message; that there is a class of states, related to SU(3) 
configurations of maximum weight, that can gain huge 
amounts of correlation energy provided the structure of 
the spherical mean field does not hamper their collective 
build-up. In the present case, instead of the assumed 
SU(3) coupling scheme in both shells, the physical situa- 
tion is closer to pseudo-SU(3) [2l| in the sd shell, based in 
the quasi-degenerate 51/2-^3/2 doublet and quasi-SU(3) 
(2^ in the pf shell, originated in the Aj=2, Al=2 pair of 
orbits, /7/2-P3/2- This means lower symmetry and there- 
fore less correlation energy. But this reduction is more 
than compensated by using a better spherical mean field, 
that, in particular, must incorporate quadratic monopole 
terms. These, as we shall see, reduce notably the mean 
field excitation energy of the np-nh configurations 23^. 
Furthermore, pairing correlations need to be taken into 
account as well. 

In addition to the estimates of the energy gains, the 
above SU(3)-like models make it possible to compute the 
intrinsic quadrupole moments associated to the differ- 
ent np-np deformed structures. Using the pseudo-|-quasi- 
SU(3) prescription we find Qo=125 e fm^ for the 4p-4h 
configuration and Qo=180 e fm'^ for the 8p-8h, consis- 
tent with what one would expect for a deformed and 
a superdeformed band. Let's underline that what we 
submit is that super-deformation in this region of nu- 
clei can be achieved readily in the space of two major 
oscillator shells, provided that SU(3)-like geometries of 
the spherical mean field are at hand. The SU(3) limit 
gives Qo=148 e fm^ and Qo=226 e fm^ instead. As 
the quadrupole correlation energy should vary as the 
square of the quadrupole moment, we can reckon that the 
pseudo+quasi-SU(3) correlation energy represents 2/3 of 
the SU(3) limit. 



III. VALENCE SPACE AND EFFECTIVE 
INTERACTION 

As we have explained above, an adequate valence space 
for the study of the coexisting bands in "'"Ca consists of 
the sd and pf major oscillator shells. Its only drawback 
is that the dimensions involved (10^^) are beyond our 
present computing capability. A possible way out is to 
close the dt^/2 orbit, that is, to work with a virtual ^^Si 
inert core. The quadrupole colectivity of the solutions 
will be reduced by this truncation, but we have checked 
in ■^^Ar that the effect is moderate. As a bonus, this 
truncation reduces drastically the spurious center of mass 
components of our wave functions that can therefore be 
controlled perturbatively. Hence, our valence space will 
encompass the lsi/2, Od3/2, 0f7/2, IP3/2, Of5/2 and lpi/2 
orbits, leading to maximum basis sizes of O(IO^). Fol- 
lowing the notation of ref. [l^] we call this valence space 
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r2pf ■ In an harmonic oscillator major shell of principal 
quantum number p, represents all the orbits except 
the one with the largest j = p + \- 

For this valence space we will use the effective interac- 
tion of ref.Q, named SDPF.SM and described in detail 
in [iJ]- In addition, to reduce the effects of the mixing 
with states with spurious centre of mass excitations, we 
add to the interaction the center of mass hamiltonian, 
^cmHcm, with XcM = 0.5. The effect of this correc- 
tion is small, since the blocking of the ^5/2 orbit reduces 
greatly the centre of mass contamination. An unwelcome 
effect of the blocking of the 1^5/2 orbit, we have to live 
with, is that it does not affect equally to the different 
np-nh configurations. It has none on the closed shell and 
becomes maximal for the 12p-12h excitations. In ref.Q 
the losses of correlation energy were compensated by an 
attractive monopole term parabolic in n. We shall do the 
same in the present case, but, instead of trying to infer 
from partial calculations the parameters of the monopole 
correction, we shall adjust them so as to locate the two 
excited 0+ states close to their experimental values. We 
discuss this issue in more detail along with the presenta- 
tion of the results of the fully mixed calculations. 




Figure 1: Superdeformed band in Ca. _E-,'s, experimental 
vs fixed 8huj results 



Table I: Quadrupole properties (in efm^ and e^fm^) of the 
superdeformed band in ''"Ca in a fixed 8p-8h calculation. 



IV. CALCULATIONS AT FIXED Nhu 

With the effective interaction and the valence space 
ready, we start by making calculations at fixed nhuj . 
The aim is to verify that our expectations based in pre- 
vious experiences and in algebraic models are fulfilled. 
And indeed they are. In fig[T] we present the results of 
the 8huj calculation compared with the experimental su- 
perdeformed band from ref. J^] in a backbending plot. 
The accord is excellent, only the slight change of slope at 
J = 10 is not reproduced by the calculation. Notice that 
the band is very regular, showing no backbending up to 
J = 24, contrary to the situation in ''^Cr or "^^Ar. This 
delay in the alignement regime is surely due to the extra 
collectivity induced by the presence of eight particles in 
quasi-S't/(3) orbitals -four in ^^Ar- and four particles in 
pseudo-S'[/(3) orbitals- absent in "'^Cr. The tiny back- 
bending at J=20 produced by the truncated calculation 
of ref. [isj . disappears when the unrestricted r2pf va- 
lence space is employed. 

In table |I] we have collected the quadrupole proper- 
ties of the superdeformed band. We have computed 
the B{E2)^s and the spectroscopic quadrupole moments 
Qspec using the conventional effective charges Sqt^ = 
6q„ = 0.5 and b=1.974 fm. We extract the intrinsic (tran- 
sition) quadrupole moments Qo{t) from the B{E2ys and 
the static ones Qo{s) from the spectroscopic quadrupole 
moments, assuming K = and using the standard for- 
mulas relating laboratory and intrinsic frame quantities; 



J 


Qspec. 


B(E2)j^j_2 


Qo(s) Qo(t) 


2 


-51.4 


643 


180 180 


4 


-64.6 


905 


178 178 


6 


-68.4 


968 


171 176 


8 


-69.6 


980 


165 173 


10 


-69.9 


953 


161 168 


12 


-70.5 


896 


159 162 


14 


-72.0 


804 


159 152 


16 


-72.5 


696 


159 141 



B{E2, J ^ J - 2) = ^\{JK2Q\J - 2, K)\'' Ql{t) (3) 
K^l/2, 1; 

The value of the intrinsic quadrupole moment roughly 
corresponds to a deformation parameter /3 = 0.6, which is 
characteristic of a superdeformed shape and agrees with 
the experimental value of ref.p^ (a comparison with a 
subsequent analysis of the same experiment 25] will be 
made in sections IVI Al and I VI Cp . In Fig. [2] we have plot- 
ted the calculated and experimental results; they agree 
within the large experimental error bars. Notice that 
as J grows, the theoretical results lose some collectiv- 
ity, whereas the experimental fit to the Doppler Shift 
Attenuation data is compatible with a constant transi- 
tion quadrupole moment. The experimental point corre- 
sponding to the 4+^2+ transition comes from an earlier 
measure of the lifetime of the 4+ state and its in-band 
branching ratio [1^ . 

The rather impressive agreement of the energetics of 
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Table II: Comparison of the experimental transition energies 
(in keV) in the 7 band of ''"Ca with the fixed 4huj calculation. 

Transition E^(4p-4h) (Experiment) 



1 
4+ 

7+ 
9+ 

13+ 



1 
7+ 



11; 



819 

1244 

1187 

1501 

2346 

1518 

1943 



781 

1260 

1369 

1538 

2773 

1827 

3044 



100 



8 10 
J 



12 14 



16 



Table III: Quadrupole properties (in efm^ and e^fm'') of the 
ND of *"Ca in a fixed Ahtu calculation. 



Figure 2: Superdeformed band in Ca; transition quadrupole 
moments Qo(t); experimental results from ref. [23| (lozenges) 
and from ref. [101 (circles) vs the results of the fixed 8p-8h 
calculation (squares). 



J 




(MeV) 

Figure 3: ND band in *'^Ca. iS^'s, experiment vs 4p-4h cal- 
culation. 



the superdeformed band with the calculation at fixed 
8hio is reminiscent of what we had found in the study 
of ^^Ar and, as it was the case then, suggests that the 
SD band is of rather pure 8p-8h character. 

We have performed as well a fixed 4p-4h calculation in 
order to describe the ND bands in ^°Ca. The most salient 
aspect of the calculated results is the triaxial character 
of the solution, with a well developped 7-band based in 
the second 2+ state (2+) of the 4p-4h configuration. The 
results for the band based on the lowest 0"*" state of the 
4fia' space are compared to the experimental ones from 
ref. [lOj (band 2) in Fig. [3) The agreement is now much 
worse that for the SD band, although the main trends 
are already present in the calculation. At J = 12 the ex- 



J 


Qspec. 


B(E2),7^j_ 


-2 Qo(s) QoW 


2 


-31.0 


269 


109 116 


4 


-41.0 


364 


113 113 


6 


-48.9 


341 


120 104 


8 


-44.0 


309 


104 97 


10 


-48.9 


237 


113 84 


12 


-38.3 


86 


86 50 


14 


-40.9 


115 


91 58 


16 


-34.8 


47 


76 37 



perimental band upbends while the calculation produces 
a strong backbending. 

The 2+ band-head of the 7-band is located 2.05 MeV 
above the 0+ band-head of the K=0 ND band. Experi- 
mentally the splitting is 1.90 MeV. The transition ener- 
gies inside the 7-band are compared to the experimental 
ones, ref. [l3| (band 4), in table [Til The calculated values 
compare reasonably well with the experiment, except for 
the location of the 13+ state that appears 1 MeV too low. 

The quadrupole properties of the K=0 ND band are 
shown in tabic IIIll In this case, the corresponding defor- 
mation parameter has the value f3 = 0.3, a typical value 
for a normal deformed band. We postpone the compar- 
ison with the experimental data to the discussion of the 
complete calculation. 

We close this section by stressing that the shell model 
calculations in the r2pf valence space at fixed 8p-8h and 
4p-4h configurations give a convincing description of the 
superdeformed band in "'"Ca, and are also able to de- 
scribe the corresponding ND rotational spectra, which 
turns out to be that of a triaxial structure developing a 
K=0 and a K=2 (7) band, in agreement with the exper- 
imental information. Indeed, in the physical states the 
configurations with different values of nhu; mix, mainly 
through the cross shell pairing interactions, and it will be 
the task of the next sections to understand how this mix- 
ing proceeds and to compare the complete results with 
the experimental data. 
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V. THE ENERGIES OF THE Nfto; 
BAND-HEADS RELATIVE TO THE CLOSED 
SHELL 



Before moving into the full r2pf space diagonaliza- 
tions, we need to explore more in depth the informa- 
tion gathered in the fixed Nftw calculations. With the 
SDPF.SM interaction, the nhui bandheads lie too high 
in energy relative to the closed shell, as it was the case 
in our study of '^^Ar. A small part of this missing energy 
can be due to residual defects of the cross-shell monopole 
terms of the interaction, but the bulk of it is due to the 
blocking of the 0d^f2 orbit. We can absorb this effect via 
the modification of two global monopole terms: 



A(?i) = -n{Ti - l)6i + n{12 - n)&2 



(4) 



where n is the number of particles in the p/-shell, and the 
(5's are the modifications of the global p/-shell monopole 
interaction and the r2-pf monopole interaction, respec- 
tively. In our study of '^^Ar, an equivalent expression 
was used. In a first set of exploratory calculations, we 
realized that when the S's were chosen so as to locate 
the two excited 0^ states at about their experimental 
excitation energies, the percentage of closed shell in the 
ground state of ''"Ca was too low (54%), i. e. the mix- 
ing was too strong. In order to diminish the mixing we 
have multiplied all the cross shell off-diagonal matrix el- 
ements by a factor 0.8. Now, the values of the J's that 
place the 0+ states correctly are 5i=-0.27 MeV, and 82=- 
0.130 MeV. With this choice the closed shell component 
of the ground state rises to 65%, that we take as a rea- 
sonable value. We have used another set of S's that gives 
75% of closed shell and no dramatic changes are observed 
in the ensemble of observables. As the 0.8 scaling implies 
a reduction of the pairing interactions in the space, the 
isovector pairing matrix elements of the p/-shell part of 
the SDPF.SM interaction, that were reduced by a factor 
0.6 in our study of '^^Ar, are taken now as 0.8 times their 
SDPF.SM values. In summary, we have loosely fitted 
two global monopoles of the interaction to the location 
of the excited 0+ states, keeping the core excited com- 
ponents in the 30% range. Actually, the difference in 
correlation energy between a full sd-shell calculation for 
the nucleus N = Z = 20 — ^, and the calculation with 
the orbit Od5/2 blocked, accounts for two thirds of the 
value of A(n). The pairing modifications can be justified 
with perturbation theory arguments (see the discussion 
in [13). 

We have now all the ingredients needed to make quan- 
titative the discussion of sec. [ITl They are plotted in 
Fig. m The squares give the lowest uncorrelated energy 
in each NHlo space (to be more precise, the minimum of 
the expectation value of the SDPF.SM Hamiltonian plus 
the correction of equation [H calculated for all the Slater 
Determinants of the m-scheme spherical basis). As we 
had anticipated the uncorrelated energies do not grow 
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Figure 4: Energies in the different np-nh configurations; Low- 
est Slater Determinant plus the correction of equation 3] 
(squares); Energy of the lowest 0"*" state in the r2pf space 
at fixed Nfioj (circles). The diamonds give the energies of the 
fully mixed calculations in the r2pf space 



linearly with the number of particles excited across the 
N=Z=20 gap (if it had been so, all the squares would lie 
in the prolongation of the straight line joining the n=0 
and n=2 points). The increase is much slower, with a su- 
perimposed odd-even effect in n/2. Notice that energy of 
the lowest Shut Slater Determinant lies « 20 MeV above 
the closed shell. Finally, we carry out the explicit diago- 
nalizations separately for each Nhio value and we find the 
results plotted as circles. We can see that the correlation 
energies are very large, in particular for the SD band- 
head that gains 18.5 MeV and becomes almost degener- 
ated with the closed shell configuration. We can make a 
rough analysis of the correlation energies in terms of their 
T=0 and T=l pairing and multipole-multipolc (mainly 
quadrupole-quadrupole) content. For that we compute 
the expectation value of the (monopole free) pairing part 
of the effective interaction in the fixed nhuj bandheads 
and subtract it from the total correlation energy. In the 
4huj 0"*" ND state, out of the 9.5 MeV of correlation en- 
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2 4 6 8 10 
np-nh configuration 

Figure 5: Energies of the yrast states in the fixed np-nh cal- 
culations (corresponding to the circles in Fig. [3] for J=0) 



ergy, 5.5 MeV come from T=l pairing, 0.5 MeV from 
T=0 pairing and 3.5 MeV from the quadrupole correla- 
tions. In the %fujj 0+ SD state, the share is 5.5 MeV, 
0.5 MeV and 12.5 MeV respectively. Thus, the contribu- 
tion of the neutron proton pairing amounts to 2.33 MeV 
for any of the bands. 

The structure of the line joining the black circles in the 
plot in quite interesting, because it resembles the energy 
versus deformation curves typical of projected Hartree- 
Fock calculations. The number of particle-hole excita- 
tions in the X-axis can be taken as roughly proportional 
to the deformation. With this perspective, we distinguish 
three coexisting minima, spherical, n=0, deformed, n=4, 
and superdeformed, n=8, separated by the n=2 and n=6 
maxima. The n=10 band-head lies very high in energy 
and we have verified that it plays no role at all in this 
problem. The structure of this "potential energy" curve, 
consisting of wells and barriers, provides an intuitive pic- 
ture of how the deformed and superdeformed bands can 
preserve their identity after mixing. 

Keeping on with the mean field analogy, it is worth 
to note that the black circles in Fig. 2] correspond to 



Table IV: Percentage of np-nh components and energy of the 
first three 0+states (GS, ND, and SD) of ""Ca 





Op-Oh 2p-2h 4p-4h 6p-6h i 


Bp-8h 


E(th) E(exp) 




65 29 5 







0+ 


1 1 64 25 


9 


3.49 3.35 




9 4 


87 


4.80 5.21 



J=0^ states, i. e. they are angular momentum pro- 
jected solutions (before variation, VAP in the usual jar- 
gon). In Fig. Owe add the points for all the lowest states 
of even J that appear in the deformed {Afku and Ghuj ) 
and superdeformed {8huj ) bands and the corresponding 
2huj and lOhio states. We can observe that the structure 
of wells and barriers of the J=0 curve, that protects the 
deformed and superdeformed bands from strong mixing, 
is maintained up to J=10. Above, the barriers flatten 
and we should expect larger mixing. In particular this 
can lead to the erosion of the collectivity of the SD band 
at high spin in the calculation. 

VI. RESULTS OF THE COMPLETE 
CALCULATIONS IN THE rap/ SPACE 

The calculations in the r2pf space are computation- 
ally very demanding. In addition to the very large di- 
mensions involved, the need to compute several states of 
the same angular momentum increases substantially the 
number of Lanczos iterations needed to achieve conver- 
gence. The shell model code ANTOINE has been used 
throughout. 

The structure of the first three O"*" states, the spher- 
ical ground state and the excited deformed and super- 
deformed band-heads is shown in table IIVI The two 
body cross-shell off diagonal matrix elements can con- 
nect directly configurations differing only in two particle 
hole jumps. We see in the table that, indeed, the ground 
state Op-Oh leading component is mainly correlated by 
2p-2h components. The leading 4p-4h component of the 
ND state can, in principle mix with 2p-2h and 6p-6h 
components. Actually the mixing is dominated by the 
latter, with a non-negligible 8p-8h piece. The leading 
8p-8h component of the SD could also mix with 6p-6h 
and lOp-lOh configurations directly, but it chooses none. 
The SD band-head is very pure, with only small amounts 
of 4p-4h and 6p-6h components. Let us mention that the 
yrast band in the 6hiu space, corresponds also to a rotor 
with deformation slightly larger than the one in the Ahuj 
space. On the contrary, the deformation and the corre- 
lation energy in the lOhiu space are smaller, leading to 
negligible mixing with the other spaces. 

In the process of mixing the winner is, energy wise, the 
ground state that gains almost 5 MeV, mostly pairing- 
like. The ND band gains barely 2 MeV and the SD band 
essentially nothing. That's why, in order to reproduce the 
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Figure 6: .E^-'s in the SD band, mixed calculation vs. experi- 
ment 

experimental situation, the three 0'*' states before mixing 
must be degenerated or even with their energies inverted 
The energies of the three physical 0+ states after mixing 
are represented by the diamonds in Fig. [H 



A. The Super-Deformed Band 

We start the discussion of the results with the excita- 
tion energies of the SD band as produced by the fully 
mixed calculation. As we can see in FigjSl the mixing 
does not modify noticeably the features already present 
in the fixed 8p-8h calculation. Perhaps one could de- 
tect some irregularities at the upper part of the band. 
The calculated sequence crosses the experimental one 
at around the rotational frequency where the calculated 
states start loosing collectivity, but in global terms, the 
agreement is excellent. The results of the mixed calcula- 
tion beyond J=16 are equivalent to that of the fixed 8hu! 
calculation shown in Fig. [T] 

In table |Vl we give the np-nh structure of the SD band 
up to J=16. The 8p-8h component is nicely constant up 
to J=10, as we expected from the analysis of the potential 
energy curves in Fig[51 At higher angular momentum, the 
mixing with 6huj components becomes stronger, and we 
should expect a less collective behavior. It is interesting 
to follow the evolution of the location of the SD states 
with increasing angular momentum. The J=0 SD state 
is the third J=0 state in the spectrum. The 3—2, 4, and 
6 states are the second of their spins, J=8 and 10, third, 
J=12, fourth, J=14 sixth, and finally 3—16 fourth again. 

We shall examine now the quadrupole properties of 
the superdeformed band in ^°Ca as they come out of the 
full r2pf calculation. In table IVII we have gathered the 
theoretical spectroscopic quadrupole moments and the 
B(E2) values. We can conclude that the mixing causes 
just an erosion of the 8huj values presented in table U up 



Table V: np-nh structure of the superdeformed band of Ca. 



J Op-Oh 2p-2h 4p-4h 6p-6h 8p-8h 






9 


4 


87 


2 


11 


4 


85 


4 


8 


5 


87 


6 


3 


5 


91 


8 


2 


6 


91 


10 - 


1 


12 


87 


12 - 


2 


29 


69 


14 - 


11 


27 


63 


16 - 





40 


60 



Table VI: Quadrupole properties (in efm^ and e^fm*) of the 
superdeformed band in ''"Ca. Full r2pf calculation. 



J 


Qspec. 


B(E2)j^j_ 


-2 Qo(s) Qo(t) 


2 


-47.7 


579 


167 171 


4 


-61.1 


813 


168 169 


6 


-66.3 


874 


166 167 


8 


-66.5 


906 


158 166 


10 


-66.3 


844 


153 158 


12 


-71.8 


546 


162 126 


14 


-62.1 


557 


139 127 


16 


-64.9 


429 


142 111 



to J=10-12. Beyond, the larger presence of less collec- 
tive 6huj components has a much stronger effect, partic- 
ularly in the B(E2)'s. Thus, the values of the transition 
quadrupole moments diminish rapidly at the end of the 
band. The static quadrupole moments vary less abruptly. 

In Fig. [7] we compare our calculated numbers with 
the experimental results. Comparing with the fixed 
8huj results, we find a large reduction of the transition 
quadrupole moments of the three upper transitions of the 
band. It is possible that this drop in collectivity be an 
artifact of the r2pf space, but for the moment we have 
no way to verify it. 

The experimental points in the figure for J>6 are those 
of the first analysis of the Doppler Shift Attenuation 
data made in ref.fTo"!. The 3=4 point comes from an- 
other, earlier, experiment 26]. The calculated transition 
quadrupole moments agree with the data for J<10, but 
underestimate them for the three uppermost transitions. 
In rcf. 25] a reanalysis of the same experimental data was 
made, and the preferred solution differs from the prece- 
dent one. Whereas the former solution was a constant 
Qo(t)=180to'29 e fm^ for the six transitions from the 
states of the SD band with J=16 to J=6, the latter gives 
Qo(t)=181+°;26±0.21 e fm^ for the J=16 to J=12 states 
and Qo(t)=118j;°;°^±0.13 e fm^ for J=10 to J=6. No- 
tice that, compared to ref. 10] an extra systematic error 
has been added. In spite of that, the final error bars of 
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Figure 7: Transition quadrupole moments in the SD band. 
Full r2pf results (squares) compared to the experimental data 
from refs. [26l ] (lozenges) and (circles) 

the values corresponding to the lower J transitions are 
largely reduced with respect to those of the upper tran- 
sitions. Actually, their extracted Qo(t) values are barely 
compatible with the value for the 4"'"-^2+ 914 keV tran- 
sition, Qo(t)=169t!]:32 e fm^, obtained in ref. ^ (the 
lower and upper tips of one and another error bars just 
touch at the value 137 e fm^). 

If the results of the analysis of Chiara et. al. were 
the only viable interpretation of the data, the compari- 
son with our calculated Qo(t)'s. would be rather poor, 
meaning that something important in missing in our ap- 
proach. The beautiful agreement we had for J<10 will be 
lost, and the fact that the increase of the experimental er- 
ror bars for the upper transitions make our results deviate 
less, is a meager compensation. We have tested solutions 
with larger mixing, through the mechanisms discussed in 
section |Vl but even going as far as keeping just 50% of 
closed shell in the ground state, we obtain only an extra 
10% erosion of the B(E2)'s for J<10. On the contrary, 
more mixing brings in large reductions of the transitions 
in the upper part of the SD band. In the section dealing 
with the out-band decay branches we propose a mech- 
anism that conciliates the theoretical picture with the 
experimental data. It is based on a detailed analysis of 
the decay of the J=8 member of the SD band, a decay 
that in our opinion, has driven the fit of ref.[25| into the 
new set of Qo(t)'s. 

B. The Normal-Deformed Triaxial Bands 

We move now to the full r2pf space results for the two 
bands that according to our calculations have a dominant 
4p-4h structure. As we have shown in section IIVI they 
are consistent with the presence of a deformed triaxial 
intrinsic state. In table IVIII we list the percentage of 



Table VII: np-nh structure of the deformed band of Ca 
based on the first excited 0^ state 



J Op-Oh 2p-2h 4p-4h 6p-6h 8p-8h 



1 


1 


64 


25 


9 


2 


1 


64 


24 


10 


4 


1 


68 


23 


8 


6 


2 


75 


20 


3 


8 


21 


62 


15 


2 


10 - 




81 


17 


1 


12 - 




81 


18 


1 


14 - 




82 


17 


1 


16 - 




79 


19 


1 



the different np-nh components in the deformed band of 
^"^Ca based on the first excited 0^ state. We observe 
that the 4p-4h dominance is less strong than the 8p-8h 
dominance in the SD band at low spins, and larger at high 
spins. Another characteristic feature is that the mixing 
proceeds through the 6p-6h components, with the 2p-2h 
components completely absent except in the J=8 state, 
where an accidental degeneracy occurs. We do not have 
a definite explanation for this preference, that can be due 
to phase space considerations (there are much more 6p-6h 
states to mix with than 2p-2h states) but more probably 
to the fact that the collectivity of the 4p-4h and 6p-6h 
spaces is very similar and much larger than that of the 
2p-2h space. That the mixing strength of the 6p-6h space 
is exhausted by the ND band could explain why the SD 
band is so pure. As expected in a collective picture, the 
spread of the wave functions of the states of the 7 band 
among the np-nh spaces is very similar to that of the ND 
band. 

In Fig. [5] we plot the calculated energies of the ND 
band compared with the experiment. We observe that 
the mixing improves clearly the agreement in the lower 
part of the band. The discrepancy at the backbending is 
the same already present in the Ahio calculation. 

The 2+ band-head of the 7 band is predicted at 
5.88 MeV of excitation energy compared with the ex- 
perimental value 5.25 MeV. This means that the mixing 
increases the splitting of the ND and 7 bands by 350 keV 
with respect to the result of the 4Sw calculation. The in- 
band excitation energies change very little with respect 
to the ^hhj values gathered in table [III 

The quadrupole properties of the ND band -table 
IVIIII - are very similar to those of the '^fiio calculation. 
Only at J=8 there is a reduction in collectivity due to 
the accidental mixing with a nearby 2p-2h state that we 
have already commented. As can be seen in Fig. [51 where 
we have plotted the transition quadrupole moments, this 
reduction goes in the direction demanded by the data. 
Globally the agreement is quite good, the trends are very 
well reproduced and in most cases the theoretical num- 
bers fall inside the experimental error bars. At the upper 
part of the band the theoretical values underestimate the 
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Figure 8: Deformed band in Ca. E-,'s, experiment vs full 
r2pf calculation. 



Table VIII: Quadrupole properties (in efm^ and e^fm*) of the 
ND ^''Ca, r2pf calculation 



J 


Qspec. 


B(E2)j^j_2 


Qo(s) Qo(t) 


2 


-32.2 


292 


113 121 


4 


-42.1 


397 


116 118 


6 


-47.3 


346 


118 105 


8 


-35.5 


227 


84 83 


10 


-48.1 


161 


111 69 


12 


-37.6 


75 


85 47 


14 


-39.1 


112 


87 57 


16 


-35.6 


49 


78 37 
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Figure 9: Transition quadrupole moments in the SD band. 
Full r2pf results (squares) compared to the experimental data 
from refs. [131 (lozenges) and [13] (circles) 



Table IX: Quadrupole properties (in efm^ and e^fm*) of the 
7 band in ^"Ca, r2pf calculation. 



J 


Qspec. 


B(E2)j^j_ 


_2 B(E2),,^,,_ 


-1 Qo(s) 


Qo(t) 


2 


28.6 






100 




3 


-0.24 




427 






4 


-17.7 


133 


284 


122 


106 


5 


-25.7 


214 




111 


106 


7 


-33.5 


211 




104 


90 


9 


-35.6 


175 




96 


77 


11 


-45.9 


110 




115 


59 


13 


-40.5 


80 




97 


49 



experimental ones. The data come from different sources, 
for the J>8 they are taken from the Doppler Shift Atten- 
uation analysis of ref. [10], for J=6 and J=8 from the 
lifetimes and branching ratios measured in refs.[2^.[27|. 

In table llXI we have collected the calculated quadrupole 
properties of the 7 band. The intrinsic information 
has been extracted from the B(E2)'s and spectroscopic 
quadrupole moments assuming K=2. We find again a 
well defined deformed intrinsic state, with similar static 
and transition quadrupole moments, both very similar to 
those of the ND band. 



C. Out-band transitions 

The experimental data on out-band transitions are 
scarce and very often affected of large uncertainties. For 
the low energy part of the spectrum of '^"Ca lifetimes and 
branching ratios are known for some levels. These will 
be compared with our predictions in tables |X] and IXIII 
For J>6, we rely in the semi-quantitative informations of 
refs. [l3] and [23. In table |Xl we focus in the out-band 



transitions from ND and SD states. The balance is un- 
even. The transition probabilities of the decays from the 
2+ states in the ND and SD bands to the ground state 
are largely under-predicted by the calculation. When we 
increase the mixing they increase at most by a factor 
two, very far from what would square with the data. For 
these small B(E2)'s it is possibly not sensible to reason 
in terms of factors, but instead think that some additive 
contribution is lacking in our space. It is worth to recall 
here that the span of B(E2) values that we are trying to 
explain simultaneously is of three orders of magnitude. 
The decay of the superdeformed 0+ is well reproduced, 
and the same applies to the decay of the superdeformed 
4+ to the deformed 2+. The in-band transition is simul- 
taneously well accounted for (see Fig. [7]), but apparently 
we miss badly the two remaining transitions known ex- 
perimentally. We cannot but recognize our limitations. 

In table IXII we present the calculated B(E2)'s, E^'s 
and the in-band branching ratios. For the experimental 
information we draw from the article by Chiara et. al. 
[2^ . The calculated branching ratios for the three up- 
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Table X: Out-band transitions from the superdeformed (SD) 
and normal deformed (ND) bands of ''"Ca. The energies are 
in keV and the B(E2)s in e^fm"*. 



Transition 



2+ 



B(E2) 



Theory Experiment Theory Experiment 



0^5 


4232 


3904 


1.8 


18±1 


2+ 


565 


1307 


58 


134.8±24.5 


0^5 


5263 


5629 


0.1 


1.7±0.4 


0+ 


1769 


2277 


3 


20.9±5.0 


2+ 


2045 


2638 


19.4 


21±4 


^ND 


882 


1264 


6.8 


116±34 


2+ 

^3 


397 


1294 


2.7 


176±41 



Table XI: Comparison between the theoretical in-band 
{SD — > SD) and out-band transition probabilities {SD 
ND) for the states of the superdeformed band of ''"Ca with 
J>6. The energies are in keV and the B(E2)s in e^fm'*. 



Transition 


B(E2)j^j_2 


E^ 


E.^(th) BR % 




> 4+ 

4d 


874 


1432 


1521 




43 


2695 


2403 


46 


^SD ^ 


' 6+^ 


906 


1880 


2015 






6+ 


7.8 


2921 


2904 


93 




^SD 


844 


2481 


2371 




S+ 
°ND 


12 


3030 


2929 


96 


12+ - 


- 10+^ 


546 


2932 


2857 






1.5 


3.590 


3750 


100 




-+ 12+ 


557 


3230 


3100 




12+ 


5.0 


4264 


3952 


97 




- 14+D 


429 


3563 


3447 




14^D 


0.03 


5531 


5520 


100 



permost transitions of the SD band are close to 100% in 
agreement with the experimental observation. For the 
next transition, the predicted 96% branchin g ra tio looks 
too large when inspecting the figure in ref. [2^ , but the 
paper does not give a figure for it. However, it is in the 
decay of the 8+ member of the superdeformed bands that 
we depart dramatically from the experimental branching 
ratio, 96% calculated vs 20% experimental. We shall de- 
vote some space to this comparison, mainly because it is 
in our opinion due to this very number that the authors 
of ref. ^25^ obtain a fit to their Doppler Shift Attenuation 
data in which, for J<10, the transition quadrupole mo- 
ments of the band have values that correspond actually 
to a normally deformed band. Before that, let's men- 
tion that the branching ratios of the 6+ SD state seem 
also consistent with the data of Chiara et al. and with 
a rather pure superdeformed character. We have already 
seen that the calculation also describes correctly the 4+ 
SD branching ratios. 



Our argument goes as follows: Given that the phase 
space factors favor the out-band transition by a factor 
9, and considering that the B(E2)'s of the 8g^ — »■ 6gjj 
and S~l^jj ^~ND tables [1] and IIIII are 980 e^fm''^ and 
309 e^fm^ respectively, it is readily seen that the experi- 
mental branching ratio cannot be reproduced unless the 
9307 keV and 9856 keV experimental states correspond 
to a 50% mixing of the pure SD and ND states. Assum- 
ing that the J=6 states are pure ND and SD, this leads 
to BR=27%, with in-band B(E2)= 490 e^fm" and out- 
band B(E2)= 150 e^fm^. The in-band transition in the 
ND band should have also B(E2)= 150 e^fm'^. If we fur- 
ther assume that the J=10, 12338 keV, state of the SD 
band is pure, it will decay equally to the 9307 keV and 
9856 keV J=8 states, which seems to be the experimen- 
tal situation. Translating this into transition quadrupole 
moments, we should have Qo(t)(10+)=118 e fm^ and 
Qo(t)(8+)=122 e fm^ in the super-deformed band and 
Qo(t)(8+)=69 e fm^ in the normally deformed band, in 
excellent agreement with the experimental analysis of 
refs. H^- As a bonus, the low spin part of the band 
remains truly superdeformed. Actually the band is su- 
perdeformed all along except for the 8+ state, even if the 
accidental degeneracy of the ND and SD 8+ states pro- 
vokes a strong reduction of the B(E2)'s of the transitions 
to and from the latter state. In addition, the lifetime of 
the 9856 keV, J=8, state is reduced by about a factor 
two. 

Then, why are these features absent in our calcula- 
tion? It is clear that in order to obtain a 50% mixing 
of the ND and SD states, they must be degenerate be- 
fore mixing, their effective splitting being not larger than 
about 200 keV. To match this requirement is beyond the 
accuracy of our theoretical description. Paradoxically, 
the calculated excitation energies of the 9307 keV and 
9856 keV J=8 states, 9260 keV and 9820 keV look aston- 
ishingly precise. But we fail to have the first excited 8+ 
state, experimentally at 8103 keV at the right energy; it 
its predicted at 8900 keV. This state has a 2p-2h aligned 
nature and the fact that it mixes strongly with the ND 
state, that we have already discussed, means that both 
are degenerate before mixing at about their mean exci- 
tation energy. 

In a sense, the character of the SD band is closer to 
what is suggested by the calculated static quadrupole 
moments in table IVII The only modification brought in 
by the 8+ anomaly would be a reduction of its Qo(s) from 
165 e fm2 to 135 e fm^ 

Finally, in table IXIII we have collected the out band 
transitions of the low spin states of the ND 7 band. The 
experimental information is taken from ref. [27]. Our 
first concern is to figure out to which extent our hypoth- 
esis of low K-mixing is correct in the calculation. For 
that we compare the 2+ — > 0^^ and the 2^^ 0^^ 
B(E2)'s, the former being twenty times smaller. A simi- 



lar reduction is found in the 4+ 



2+ 



VS 4 



ND 



case, thus confirming the validity of our assumption. The 
7 energies compare quite well with the experimental re- 
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Table XII: Out-band transitions from the 7 band of Ca. 
The energies are in keV and the B(E2)s in e^fm*. 



Transition B(E2) 

Theory Experiment Theory Experiment 



0^5 


0.20 


1.0±0.3 


5690 


5249 


0+ 


16 


100±3 


2277 


1896 


2+ 


91 


179±50 


1575 


1343 


2+ 


28 


27±5 


2295 


2125 


2+ 


20 


32±7 


2738 


2603 


4+ 


49 


49±34 


1632 


1229 



ture. Their energetics is not very satisfactory. The ND 
band-head, that is experimentally the first excited state, 
is calculated about 5 MeV too high, while the lowest 0+ 
pertains to an oblate side band not present in the experi- 
mental data. A bunch of states of unknown nature (J=2, 
4, 6, and 8) are even lower. The SD band is also high 
but lower than the ND band. The splittings in the SD 
band look healthy, but the ND band does not resembles 
at all to a band. The 7 band appears also at very high 
energy. Concerning the transition probabilities they are 
reasonable for the ND band and rather small for the SD 
band. No out-band transitions are reported. 



suits. The B(E2)'s are in excellent agreement with the 
experiment in all but one case, the corresponding to the 
2+ — > 0^^ transition that is severely under-predicted. 
Even so, the experimental value is three time smaller 
than the equivalent transition inside the ND band. 

D. Other approaches 

There has been two other recent attempts to explain 
theoretically the structure of ''°Ca such as we know it af- 
ter the experimental findings of Ideg uchi et al. [lol | and 
Chiara et al. . Bender et al. [l3 | approach the prob- 
lem in the framework of the constrained Hartree-Fock- 
BCS approximation. This is a calculation in the intrinsic 
frame, that is perfected by means of angular momen- 
tum and particle number projections. In a further step, 
the angular momentum projected solutions correspond- 
ing to different values of the deformation are mixed us- 
ing the Generator Coordinate Method. Modern Skyrme 
parametrizations are used as effective interactions. They 
produce a ND and a SD band with a span much larger 
than the experimental one. The calculations are carried 
out only up to J=6 for the ND and J=10 for the SD band. 
They assume an axially symmetric basis, therefore pre- 
cluding the description of the experimental 7 band. Their 
ground state has a structure very similar to ours. For the 
SD band their predicted transition quadrupole moments 
are quite large, but the agreement with the data is incom- 
plete. No calculation of out-band transitions is reported. 
There is a discussion of the np-nh structure of the states 
that we find difficult to translate into our language. They 
claim that the ND and SD bands do not have a well de- 
fined np-nh structure contrary to what our results show. 
However, our percentages of np-nh components cannot be 
compared with theirs directly. The second one applies the 
techniques of Antisymmetrized Molecular Dynamics . 
They use an schematic interaction of Gaussian type. The 
calculation produces many excited bands, and, in partic- 
ular, the three bands that we have studied in this paper. 
As we do, they obtain a superdeformed band of 8p-8h na- 



VII. CONCLUSIONS 

In this article, we have shown that large scale shell 
model calculations in a valence space comprising two 
major oscillator shells, can describe the coexistence of 
a doubly magic ground state and deformed and superde- 
formed bands in ^°Ca. We have explored the role of the 
spherical mean field and the correlations in bringing the 
n-particle n-hole configurations at very low excitation en- 
ergy. We have examined the mechanisms of mixing be- 
tween different np-nh configurations, and have analyzed 
the structure of the physical states in terms of these con- 
figurations. We find that the normally deformed bands 
are dominated by 4p-4h configurations with important 
mixing of 6p-6h components. The superdeformed band is 
clearly dominated by the 8p-8h configurations with small 
6p-6h and 4p-4h contributions. We submit that for J=8 
the ND and SD bands are accidentally degenerated and 
mix at 50%. The ND intrinsic structure is triaxial and 
produces a 7 band. In the ground state, the doubly magic 
configuration amounts to 65%, mixed mainly with 2p- 
2h states. The comparison of the theoretical predictions 
with the experimental data for the energetics of the three 
above mentioned bands is excellent. We have computed 
also the electromagnetic transition probabilities for in- 
band and out-band transitions. With some exceptions, 
the agreement is also very good. In addition to the chal- 
lenge of accounting for super-deformation in a spherical 
shell model basis, i. e. in the laboratory frame, our cal- 
culations and modeling support the experimental claims 
on the existence of superdeformed bands at low excita- 
tion energy in this mass region, even in doubly magic 
40Ca. 
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